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Raman spectra of ice VII and X were investigated using strong plasma Shockwave generated by laser- induced 
breakdown (LIB) in liquid water. Simultaneously, the occurrence of the hydrogen emission lines of 656 nm 
(Ha), 486 nm (HP), 434 nm (Hy) and 410 nm (H8) was observed. At 5 X 10'^ W/cm^ optical power density, 
the O-H symmetric stretching, translational and librational modes of ice VII and a single peak at 785 cm"^ 
appeared in the spectra. The band was assigned to the Raman-active O-O mode of the monomolecular 
phase, which was the symmetric hydrogen bond of cuprite ice X. The spectra indicated that ice VII and X 
structure were formed, as the trajectory of the strong plasma Shockwave passes through the stable 
Pressure-Temperature range of ice VII and X. The Shockwave temperature and pressure were calculated by 
the Griineisen model. 

The behavior of water at conditions of extreme pressure and temperature is of profound importance to 
planetary science, geoscience, and fundamental chemistry \ Of particular interest has been the pressure- 
induced transformation to nonmolecular forms, a transition connected with hydrogen-bond symmetriza- 
tion^. Ice X, the symmetric ice structure, was detected at 44 GPa in a diamond-anvil cell by Brillouin scattering^ 
Infrared (IR) measurements have provided evidence for the transition from ice VII (or ice VIII at lower tem- 
perature^) to a phase with symmetric hydrogen bonds-ice X beginning at around 60 GPa and stable to at least 
210 GPa^"^. The melting curve showed a discontinuous change in slope at about 43 (47) GPa and 1600 (1000) K 
using Raman spectra, indicating a phase transition most likely from the ice VII to the ice X phase^'^. 

Until recently, the experimental study of water phase transition at extreme conditions was limited mainly to use 
a laser heated the static diamond anvil cell. Holmes et al. obtained a spontaneous Raman scattering from shocked 
water, the strong Shockwaves were generated by planar impactors accelerated to 2.1-3.2 km/s by a two-stage 
light-gas gun^°.We performed the experiments using a high-intensity laser-induced shockwave compression to 
obtain high-temperature and dynamic high-pressure^ LIB refers to plasma production carried out by focusing 
an intense pulsed laser beam in substances. When the laser pulse duration is longer than a few nanoseconds, 
plasma cloud evolves by directly absorbing the laser power. The rapid expansion in the ionization region creates 
an intense shock wave coming out from the plasma cloud and propagating in the medium at supersonic speed. 
The plasma region is composed of a central low refractive index distribution with n < 1 and a spherical high- 
density distribution whose distance is increasing from the origin of the shockwave expansion^^. Currently, the 
long-range order ice VIII (H2O) and the short-range order ice VII (D2O) have been obtained by LIB^^'^^. 

In this letter, an experimental study has been carried out on the dynamical process taking place in the strong 
plasma shockwave generated by an intense 12 ns (FWHM) pulse beam focused into liquid water. The occurrence 
of the hydrogen emission lines of 656 nm (Ha), 486 nm (Hj5), 434 nm (Hy) and 410 nm (HS) was observed. At 
the same time, the Raman peaks of O-H stretching vibrational region, three low-frequency Raman Stokes peaks 
and their intense anti-Stokes peaks were obtained, which present: (1) The characteristic Raman spectrum that 
developed from the symmetric hydrogen bonds of ice X at 5 X 10^^ W/cm^ optical power density; (2) One 
internal-mode Raman spectrum that belongs to the O-H stretching vibration, one translational mode T, and one 
librational mode R of ice VII. These phenomena are attributed to the ice VII and X formation in liquid water 
induced by strong plasma shockwave. 

Results 

We attempted to determine the phase transition trajectory of ice by laser-induced strong plasma shockwave. We 
found that this method yields rich in atom and molecule spectra compared to the hydrostatic pressure. 
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Representative Raman spectra in liquid water under the different 
optical power density are shown in Fig. 1. The high-frequency 0-H 
stretching vibrations in H2O ice VII form a doublet, which has been 
assigned to Vi(Ai^) and v^iEg) modes in order of increasing frequency 
(ice VII, point group of O^)^^. These modes are softening with pres- 
sures^, the softening effect can be seen from Fig. l(a, c), whose posi- 
tions are corresponding with that of H2O at the around 17 GPa and 
33 Gpa pressure^^ respectively. The frequencies of all components of 
the 0-H band increase with temperature indicating that the hydro- 
gen bond in ice VII weakens as the melting transition is approached 
and molecules acquire larger translational disorder^. The hydrogen 
emission line of Ha-656 nm was also observed (Fig. 1 (a, c)). 

Measurements in the low-frequency range reveal several import- 
ant spectrum characteristics (Fig. 1 (b)). We find two special peaks in 
the spectrum at the same optical power density, a Stokes Raman line 
and its anti-Stokes Raman line, whose Raman shift are ±785 cm"^ 
that characterizes the ice X cuprite structure^^, the mode display an 
increase in frequency as the pressure increased^^, which correspond- 
ing with the pressure should be at around 65 GPa^. A likely assign- 
ment in order of increasing frequency is T^iAig) + T^yiEg) and 
TziBig) + T^yiEg) that are unresolved bands for the translational 



bands at the around 370 cm~\ Rxy(Eg) + Rz(E2g) for the librational 
bands at around 605 cm~\ the Raman shift corresponding with 
pressure should be approximate 17 GPa^^. In Fig. 1 (d), these two 
Raman shifts are around 360 and 595 cm"^ respectively, which dis- 
play an increase in frequency as the pressure increased. Another 
spectrum characteristic is that the intensity of anti-Stokes peak is 
similar to that of Stokes peak. 

Discussion 

The Carnegie group succeeded in determining the VII -X transforma- 
tion zone in 1999, this transformation was characterized by the dis- 
appearance of the Tz(Aig) + Txy(Eg) and T^iBig) + T^yiEg) (where 
the two peaks were not resolved) lattice modes of ice VII, and pro- 
gressive occurrence of the expected Raman-active F2g mode^^. The 
lower pressure limit of the transformation was found close to the 
isobar 62 GPa^^ 

In this experiment, the Raman spectra of ice VII and X can be 
obtained, because high-intensity LIB generate the strong plasma 
Shockwave, which make a dynamic high-pressure and high-temper- 
ature field around the focal point. The mechanism of ice formation 
using LIB has been discussed by our group^^ A 1 X 10^^ W/cm^ 
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Figure 1 | Representative Raman spectra of ice VII and ice X by LIB. (a) The spectrum of O-H stretch vibration (—3080 cm~^) in ice VII and Hoc 
(656 nm) line, (b) The unresolved + Txy(Tz(Aig) + T^yiEg) and Tz(Big) + Txy(Eg)) bands around 370 cm"^ and Rxy(Eg) + Rz(B2g) bands around 
605 cm~^ belong to the translational and librational modes of the ice VII respectively. The T{F2g) peak, which characterizes the ice X cuprite structure, is 
around 785 cm 4n ice X. (c) The spectrum of O-H stretch vibration (—3070 cm~^) in ice VII and Ha (656 nm) line, (d) The unresolved + Txy(Tz(Aig) 
+ Txy{Eg) and T^iBig) + T^yiEg)) bands around 360 cm~^ and Rxy{Eg) + Rz{B2g) bands around 595 cm~^ belong to the translational and librational 
modes of the ice VII respectively. The rectangle centered at zero frequency mask peculiarity associated with the subtraction of elastically scattered 
radiation. The optical power density of Fig. 1 (a, b) is 5 X 10^^ W/cm^, and that of Fig. 1 (c, d) is 2 X 10^^ W/cm^. 
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Figure 2 | Dynamic pressure-temperature phase diagram of ice in the region VII-X. Blue closes squares, dynamic pressure at different pumping energy 
for H20y corresponding blue solid lines. Horizontal doted lines at ~2 GPa and —62 Gpa, expected liquid-VII and VII-X transition lines for H2O 
respectively. Insert: A shockwave front moving through a medium, the shockwave pulse falls more slowly than it rises A magenta horizontal and a 
magenta vertical dot lines are corresponding with the dynamic pressure and temperature of one shock rear and front at pumping energy around 500 mJ, 
which corresponds to Fig. 1 (a, b). A yellow horizontal and a yellow vertical dot lines are corresponding with the dynamic pressure and temperature of the 
yellow closed circle 2 at pumping energy around 275 mJ, which corresponds to Fig. l(c, d). These closed circles are from the experimental data. 



optical power density is enough to induce breakdown of pure water 
and some organic solvents^°'^\ First, the rising edge of the laser pulse 
whose duration is about 12 ns interacts with water, and LIB is gen- 
erated when the pulse peak arrives at the focal point. Simultaneously, 
the plasma bubble that has a diameter of hundreds of microns is also 
obtained^^'^^; Second, the plasma bubble begins to collapse and pro- 
duces strong shockwave that induce water compression, the shock- 
wave speed and the particle velocity is approximate 8.4 and 2.4 km/s, 
respectively^°'^^. The strong plasma shockwave front rises faster in 
water, whose duration is around 10 ns, and its pulse falls whose 
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Figure 3 | Experimentally determined shockwave velocity Vg and particle 
velocity Vp behind the shockwave front plotted as a function of the time 
delay between laser pulse and illumination pulse. Curve 1 and 3 show the 
shockwave velocity and the particle velocity at pumping energy around 
500 mJ, respectively. Curve 2 and 4 show that at pumping energy around 
275 mJ, respectively. 



duration is around 100 ns more slowly than it rises. Ice X can be 
formed when the shockwave front rises, the temperature is around 
900 K and the pressure is around 65 Gpa^^'^^'^^, which corresponding 
with the experimental result point 1 in Fig. 2. The water pressure 
declines slowly and uniformly^^. The shockwave generates approxi- 
mately 17 Gpa pressure and 350 K temperature^^"^" when the shock 
wave pulse falls (Point 1* in Fig. 2). Under these conditions, ice VII 
can be obtained^^. We estimated temperature by analyzing the rela- 
tive intensities of the Stokes and anti- Stokes peaks according to the 



equation T - 



hcDk 



ln-+4ln 



The intensities of 



Stokes and anti -Stokes lines depend on the population of the particles 
of the correlated ground and excited state of the specific vibrations. 
The Stokes-to-anti-Stokes intensity technique was adopted because it 
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Figure 4 | Shockwave pressures of the numerical calculations for 
pumping energy with 500 mJ (curve 1) and 275 mJ (curve 2). 
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Figure 5 | Shockwave temperature data plotted as a function of pressure. Points a and b taken to Q = 7. li^, and points c, d, e, f and g taken to Q = 8.7i?. 



is based on an intrinsic property of the studied system. The Raman 
peaks of ice X show a higher temperature than that of ice VII based on 
the equation. Third, the descending edge whose duration is approx- 
imate 12 ns of the laser pulse excites ice VII and X, and the Raman 
spectra of ice VII and X are observed. Two reasons affirm that the 
Raman spectra belong to ice VII and X. (1) The lattice translation of 
370 cm"^ at low frequency is the characteristic peak of ice VII, and 
the disappearance of the translational peak 370 cm"^ of ice VII is 
interpreted as a transition into ice X. (2) the Raman-active mode 
of 785 cm"^ characterizes the ice X cuprite structure, only the sym- 
metric hydrogen bonds can generate this Raman band in ice^^. 

The Shockwave pressure and temperature were calculated as the 
following. 

According to the plasma particle velocity Vp{t) that was obtained 
by fluorescence imaging techniques, and the relationship between Vp 

and v,^^, 



, = ci(10' 



|(Vs-Co)/c2 _ 



1) 



(1) 



V5 is the Shockwave velocity, the curves of shockwave and particle 
velocity were drew in Fig. 3. Base on the values, the shockwave 
pressure ps(^) can be also calculated using the equation, 

p, = CiPoV,(10<'''-'^°>/'^- 



-l)+pc, 



(2) 



ps{t) curves are drew in Fig. 4. According to it, the shockwave front 
rises faster, and its pulse falls more slowly than it rises. 

We put to use the Griineisen model to calculate the shockwave 
temperature^^. 



(3) 



If at a given volume, the Hugoniot pressure is known as well as the 
state (Pp T^), then Eq.(3) allows evaluation of the internal energy 
difference between these states. In this case. 
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When, at any given V, the appropriate value of y is assumed inde- 
pendent of temperature, then Py = Ps{V) can be calculated numer- 
ically, which in turn leads to the solution of Eq.(4) for the theoretical 
shockwave temperature T^. In the present case, temperature calcula- 



tions using a constant y = 0.9^^. The specific heat Q has a much 
greater effect in shockwave temperature calculations. The value of 
Cy = 7.1i^ and Cy = 8.7i^ per mol of H2O was used for low and high 
pressure, respectively. R is ideal gas constant. 

Fig. 5 shows the measured shockwave temperatures plotted as a 
function of principal Hugoniot pressure. The temperature is higher 
than that previously calculated by Stokes -to -anti- Stokes intensity 
technique, because of two reasons. (1) Water temperature drops 
quickly^^ when there are a large number of water molecules sur- 
rounding shockwave front at room temperature. (2) Shockwave 
can lead water volume V to largening rather than a constant, so 
the Hugoniot temperature is higher than actual temperature. 

In conclusion, we have observed the hydrogen emission spectra 
and Raman spectra of ice VII and X when LIB is generated in Uquid 
water. The LIB generates a high-temperature and dynamic high- 
pressure field, which lead to the intensity of anti-Stokes be close to 
that of Stokes. Ice VII and X were formed by the strong plasma 
shockwave rear and front respectively. The shockwave pressure 
and temperature are calculated by the Griineisen model. The study 
confirmed that ice X transforms from ice VII under the dynamic 
pressure condition. Last but not least we have developed a powerful 
tool to study on the phase transition of liquid system. 

Methods 

Raman scattering experiments of water. A frequency- doubled Nd:YAG laser is used 
as the excitation laser for laser- induced breakdown of liquid water, which emitted 
laser pulse duration of 12 ns with a wavelength of 532 nm. It has a repetition rate of 
1 Hz and pulse energy in the range of 0-700 mj. The laser is coupled to a 10 cm-long 
quartz cell with 25 times magnification object lens. The output light is focused on the 
AvaSpec-2048-5-USB2 fiber optic spectrometers, whose spectral resolution is 
0.05 nm, and the signals are managed by a computer. The experimental setup is 
shown elsewhere^'*. 

Fluorescence imaging experiments. The plasma fluorescence imaging of water is 
obtained by the PI- Max ICCD, whose image resolution is 1024 X 1024, the shutter 
duration is 2 ns. The interval time of sampling is 2 ns. A 656 nm filter is used to 
achieving the around 656 nm wavelength fluorescence imaging for obtaining the 
particle velocity. 

Computation. In our studies, calculations of shockwave pressure and temperature 
are performed by using the Mathematica soft. 
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